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Abstract 

The ratio between the photon spectrum in i? ^ and the differential semileptonic rate 
wrt the hadronic variable Mx/Ex is a short-distance quantity calculable in perturbation 
theory and independent of the Fermi motion of the h quark in the B meson. We present a 
NLO analysis of this ratio and show how it can be used to determine iKib/Kfel independently 
of any model for the shape function. We also discuss how this relation can be used to test 
the validity of the shape-function theory on the data. 
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1 Introduction 



The accurate determination of the parameters entering the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix V is one of the main physics goals at beauty factories. The large number 
of produced B mesons allows for the study of many decay modes and the measurement of 
branching ratios (BRs) as small as 10^^. Unfortunately, a precise extraction of the CKM 
parameters from the data is often hindered by uncertainties on the hadronic parameters 
accounting for low-energy strong interaction effects. There exist, however, observable quan- 
tities that are free, at least to some extent, from hadronic uncertainties. They are obviously 
highly valued and have been intensively looked for in recent years. 

In particular, it has been suggested many years ago that one can build a suitable ratio of 
the lepton-energy spectrum in the semileptonic decay B —>■ XJ.^ and the photon spectrum 
in B ^ Xg'^ from which the dominant non-perturbative QCD effects due to the heavy- 
quark Fermi motion drop out 0,0]. This ratio can be used to extract the combination of 
CKM matrix elements |Vub/(PtbV^*)|. Theoretically, the cancellation of long-distance effects 
takes place for hadronic invariant masses as small as Mx ~ QAqcd, where Mx is the 
invariant mass of the hadronic system and Q is the hard scale of the process (in the B rest 
frame Q = 2Ex, with Ex denoting the hadronic energy). In this regime, the dominant 
non-perturbative effects are described by a single universal function, the shape function, 
entering the different spectra. At smaller ~ ^qcd; however, few resonances dominate 
the spectra and universality is lost. It remains to be seen whether the region where the 
shape function provides a good description of hadronic physics is large enough to allow a 
clean extraction of the CKM parameters. 

Several short-distance quantities based on various partially integrated spectra or mo- 
ments have been proposed in the literature, with a particular attention to the problem of 
defining the shape function consistently beyond the leading order Alternatively, one can 
use a model for the shape function, and estimate the theoretical uncertainty by changing 
the parameters of the model. This is routinely done in order to implement cuts on Mx 
useful to identify the B X„/z/ sample 0], an approach which has been actually used 
to extract Vub from the semileptonic decay rate at LEP ||^. Very recently CLEO has also 
published a V^b measurement based on a combined analysis of the electron and photon 
spectra in semileptonic and radiative decays 0. Other proposals to identify B XJu 
events with low theoretical uncertainty include a cut on the invariant leptonic mass [[7[ 
and a combination of cuts on and Mx 0. 

In this paper, we propose a phenomenological strategy to investigate the "shape-function 
region" and possibly to extract the ratio |Kife/(Vf;,V^*)| ~ |Ktb/Kfe|- We use the fact that the 
differential rates d Trd/ dx for the radiative decay B Xg'j and d Tgi/dC, for the semileptonic 
decay B Xulv probe exactly the same long-distance structure function. The kinematical 
variables involved are x = 2E^/mB and ^ = 2t/(l + t) with t = ^1 - Ml/E\. The 
ratio of the above two distributions is a short-distance quantity, which can be computed in 
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perturbation theory. As we will explain in section 4, to a good approximation it provides a 
determination of in the SM. Including perturbative 0{as) corrections, one finds 



ub 



cb 



das)-, — (1) 



where the coefficient C{as) and the function B{S,) will be calculated in the following. 
Clearly, eq. (|I]) is valid up to higher twist 0(1 — x) ~ 0{Aqcd/ Ex) contributions. 

Our treatment of the semileptonic and radiative decays follows closely refs. [^-11]. In 



this approach, the hard scale Q is identified as twice the hadronic energy 2Ex in the B rest 
frame rather than as the heavy-quark mass mi,. Consequently, the simplest kinematical 
variable for factorizing long-distance effects in this decay turns out to be a function of 
Mx/Ex and the relevant spectrum is dTsi/dC, rather than dVgi/dMx- 
The advantages of our method are the following: 

• the cancellation of Landau-pole effects and of non-perturbative Fermi-motion effects 
in the ratio of spectra in eq. (|I]) is exact and requires no approximation; 

• the shape function is defined in a consistent subtraction scheme for infrared (IR) 
singularities, avoiding double counting of IR logs; 

• at least in principle, our ratio is valid point by point on the spectra, namely for any 
value of X = ^. 

In a realistic case, some smearing, such as partial integration over the spectra, may be 
necessary as all these calculations, including ours, rely on local quark-hadron duality. One 
could integrate eq. (0) over different kinematical ranges and check the stability of the results 
with respect to the smearing procedure. This, in turn, would amount to a test of the 
underlying assumptions, and most notably of quark-hadron duality. 

In this paper, we present a NLO calculation of the perturbative corrections to eq. (|1]). 
To make contact with experiments, we also give expressions of the relevant functions in the 
presence of kinematical cuts on the hadronic invariant mass Mx- After inclusion of the 
NLO corrections, the ratio in eq. (|l|) allows for a determination of with a 0(5%) 

theoretical error. 

The paper is organized as follows. In section ^ we introduce the formalism used through- 
out the paper, while results for the relevant spectra are given in section ^ and in the Ap- 
pendix. The extraction of |\4b/V^f,| and the test of local duality are discussed in section 
Section ^ contains our summary. 
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2 Formalism 



The heavy-quark Fermi motion in semi-inclusive decays |12|, |1^ can be treated in field 
theory by introducing the shape function /(A;+), also known as the structure or light-cone 
distribution function of heavy fiavours ||l|,Pip!^,p!5|. 

The calculation of the spectra in these decays requires the Operator Product Expansion 
(OPE) of the time-ordered product T of two transition operators (currents or effective 
Hamiltonians) . Schematically, this involves the expansion of the propagator 



Ml + 2Exk+ + k^ + ie 



Ml+ie 



+ 



Ml 



(2) 



where k+ = px ■ k/ Ex is the virtuality of the heavy quark {pi, = uibv + k) and px, Mx and 
Ex are momentum, mass and energy of the final hadronic state, including the spectator 
quarkQ. Since it is expected from the confinement dynamics that /c+ ~ ^qcdi one can 
identify three different kinematical regions, depending on the value of M^: 



1. When Ml ~ 



El S> Aq(^^, the OPE converges. At leading twist accuracy, one can 
safely retain only the lowest order (n = 0) term in eq. (@). 



Ex^QCDi all the terms in the sum become of order 1. This is precisely 



When Ml ~ 

the shape-function region. Indeed, the leading-twist terms, i.e. those obtained by 
neglecting corrections of Oik"^ /Ml) ~ 0{Aqcd/Ex), are resummed by introducing 
the shape function, formally defined as 



f{k, 



{B\hJ{k+ - iD+)K\B) 
{B\KK\B) 



(3) 



where Dj^ is the light-cone plus-component of the covariant derivative and hy is the 
heavy quark field in the Heavy Quark Effective Theory; f{k-\.) gives the probability 
of finding a heavy-quark with a light-cone residual momentum k^ inside the meson. 
Note that neglecting the terms of order k"^ /Ml in eq. (^ is equivalent to assuming 
the local quark-hadron duality. As the rate is proportional to the imaginary part of 
eq. d^), it follows that the distribution in A;+ coincides with that in —Ml/2Ex and is 
therefore experimentally accessible. 



Aqcd5 even the resummation fails and the contribution of single reso- 



When Ml 

nances to the spectra is not negligible. In this region, terms of 0{Aqcd/ Ex) are not 
negligible and the shape function cannot be defined. 



^At the leading- twist, the use of the meson mass ms instead of the quark mass mi, in the decomposition 
of heavy quark momentum pb only shifts the integration range of fc+ |lq] . 
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A different framework to factorize long-distance effects is provided by the resummation 
theory in perturbative QCD [|l^]. The connection between the two approaches is discussed 



in refs. 10|. In this framework, a factorized formula for a given distribution, for example 
the photon spectrum dVrd/dx in 5 — > Xs7, reads 

\VtbV:f [Krd{as)f{x- as) + Drd{x- a,)] , (4) 



d^rd _ T^*|2 -pO 

dx 



where Vtb and Vts are the relevant CKM elements, is a normalization factor, x is a prop- 
erly chosen factorization variable, x = 2E^/mB = 1 — M'j^/rn^^ in this example, while the 
coefficient Kj-d and the remainder function Drd{x) are short-distance quantities computable 
in perturbation theory. The form factor f{x) contains all the IR logarithms of the form 

a,(2Ex)'= (^^^yz^) "^'^^ 0<n<2k-l, (5) 

where the plus-distribution is defined as 

(6) 

1 — X J ^ 1 — X Jq ^ ~ y 

The function /(x) is well defined for 1 — x ^ AgcD/fnB- The universality of /(x) stems 
from the general properties of the soft and the coUinear radiation which produce the log- 
enhanced terms in eq. (|^). On the contrary, when 1 — x ^ KqcD/friB-, f{x) is no longer 
calculable in perturbation theory. However, /(x) is related to the shape function f{k+) in 
the following way 

fik+) = [ dk'^Cik^ - k'^; as)~f{k'^) , (7) 



where /c+ = — m5(l — x) (in the case of -B — > ^^7, the hard scale is 2Ex = fnB)- The 
coefficient function C{kj^ — k'_^;as) remains perturbative also for 1 — x ~ ^qcd/^b and 
allows for a definition of a non-perturbative form factor in this region. 



3 Results 

In this section we present a NLO determination of the short-distance functions entering the 
formulae of the photon energy spectrum in i? ^ Xg'-y and of the 2-spectrum in 5 — > XJ,v, 
written in the factorized form of eq. (^). In the latter case, we also give expressions for the 
same functions in the presence of kinematical cuts. 
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3.1 Photon spectrum in B ^ Xg-j 



The photon spectrum in the radiative B X<j7 decay can be expressed as in eq. 

S7 rate is defined as 



normahzation factor T^^ for the b 



^ rd — 



IT 



327r3 



The 



(8) 



where a = 1/137.036 is the fine structure constant. The NLO relation of the MS mass, 



"^bATsi^b) = 4.26 ± 0.09 GeV |18], with the pole mass, mj, ~ 4.7 GeV, both appearing in 
eq. (I), is 

i^b,MS i'^b) as {nib) Cf , ^( 2\ 

= 1 \- ^ ("sj . 

nib TT 

where Cp = {N'^ — 1)/(2A^) and is the number of colours. 
The coefficient function K^d reads 



(9) 



Krd{as)=Cf^ {lib) 



1 + 



+ 



E?W^(Rer.+7<?'ln^l (10) 



{lib) Cj^^ (/ib 



l^b 



2- C^Hlib) 



the definitions of rj are collected in the Appendix. We have denoted by and Cf'^ 
the LO and NLO contributions to the effective Wilson coefficients of the /S.B = 1 effective 
Hamiltonian 119,20|. The scale fib is a renormalization scale of 0{nib) and 7^7'' are elements 
of the LO anomalous dimension matrix. The terms proportional to \n{mb/fib) make the 
coefficient K^d renormalization scale-independent at the NLO. K^d depends on the top 
mass: for Mt = 174.3 ± 5.1 GeV we find 



Krd = 0.1134 [1 ± 0.014 ± 0.053;{J;|J^y 



:iii 



where the first two errors come from the top uncertainty and from the treatment of the 



charm quark mass in the matrix elements containing charm loops. We have followed |21 



and used nic/nib = 0.22 ± 0.04 in (see Appendix). The central value of eq. ( [TT| ) has been 
obtained for /if, = nib and the last error by varying fib between nib/2 and 2nib. We have also 
included in Krd a power correction — A2 {€2^^ —C[^'^ / Q) Cf^ /9m1 related to long-distance con- 
tributions from (cc) intermediate states P2[ and the electroweak effects following Ref. 



We have adopted a light Higgs boson mass, nih-, as suggested by global fits, but the sensi- 
tivity to m/j is very small. The latter two contributions modify K^d by about +2.4% and 
-3.8%. 

As is well known, NLO corrections have an important numerical impact in radiative 
decays: they increase the inclusive branching ratio by up to more than 30% In 
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eq. (|n]) the effect of NLO corrections is about (+20 ± 30)%, according to the value of fib 
between mh/2 and 2mf,. Recently, it has been shown ||2ll| that the dominant reason for the 
enhancement of the inclusive BR is the running of the b quark mass, which affects differently 
the light quarks and the top loops. Following this observation, it is possible to stabilize the 
perturbative expansion for B —>■ Xgj by a rearrangement of higher-order effects. In this 
approach, NLO corrections to Krd are not smaller, but have a different sign, 

Krd = 0.1448 (1 - 0.99 a, - 0.013) = 0.1119 [l ± 0.014 ± 0.053l[5:[!^^] . (12) 

Here the second correction is due to electroweak and 1/m^ effects and the errors have 
the same meaning as in the previous equation. Notice that now the NLO correction is 
predominantly due to the matrix elements of the effective operators. For this reason, it 
is numerically close to the NLO corrections to the coefficient in the semileptonic decay, 
as we will see later on. In particular, the large scale dependence observed in eq. (|12D is 
mostly due to the large 0{as) corrections to the matrix elements and is dominated by the 
scale dependence of ^^(/^b) and not by the scale dependence of the Wilson coefficients, in 
contrast to the situation in the preceding equation. In the following numerical studies, we 
will employ eq. (0). 

We stress that the definition of the coefficient function in eq. (0) and of the remainder 
function given below corresponds to a particular factorization scheme. Specifically, following 
the convention of |11], we absorb in f{x) only the IR logarithms of eq. (^). 



The remainder function is given by (the functions dij{x) are defined in the Appendix) 

Drd (x; a,) = Yl ^^'^ (Z^^) (Z^") ^^^^ + ^ • (^3) 

Let us comment on the behaviour of the remainder function at the end points x = 
and 1. The functions dj-/ {x) and d^s [x) have a logarithmic singularity ln(l — x) for x — >■ 1, 
coming from power-suppressed long-distance contributions. They are not factorized into 
the function f{x), which contains only terms of the form ln"(l — x)/(l — x). At the lower 
end-point x = 0, the functions dji (x) and (iyg (x) are finite, while dgg (x) has a power 
singularity dss (x — ^ 0) ~ 1/x. This limit corresponds to a soft (real) photon and the 
singularity would be regulated by including virtual photon contributions so that 1/x ^ 
(l/x) + . These contributions are not relevant to our analysis. Finally, the function dss (x) 
has a logarithmic singularity for because a non-zero photon energy regulates the 

soft but not the collinear singularity. In this case, the strange mass has to be retained (we 
adopt rris = m^/SO). The function Drd{x) is shown in fig. 1 for different values of fif,. 

We also present results for the partially integrated rate, for which a factorized formula 
also holds, 

/I 
-^dx' = \VtbV:,\^ r°, [Krd M F {x- a,) + Erd (x; a,)] . (14) 
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Figure 1: The remainder functions Drd{x) and Erd{x) for Mj 
of the renormalization scale /ij,. 



174.3 GeV and different values 



The partially integrated form factor is defined as 



F (x: a.. 



dx'f {x'] as) 



(15) 



while the cumulative remainder function is given by 



Erd (x; as) = 



dx'Drdix' 



as (yUfc 



1,8 



TT 



Note that, contrary to the dij case, the functions fij vanish at the end-point x = 1. The 
explicit expressions of the functions fij{x) are given in the Appendix, and the cumulative 
remainder function is shown in fig. 1. Considering that KrdF{x) ~ 0.1 for small x, the plot 
on the left side of fig. 1 shows that the scale dependence of the remainder function can be 
numerically important. However, the scale dependence of Erd{x) is almost entirely due to 
the coupling constant Given the expression of the partially integrated rate in eq. (0), 

we observe a significant cancellation between the 0{as) correction in eq. (p!2D and Erd{x). 
It follows, in agreement with that the NLO result for total rate is perturbatively quite 
stable. 



3.2 ^ distribution in semileptonic decay 

The starting point for dealing with spectra in B ^ Xulu is the triple-differential distribu- 



tion, which has been computed at order in ref. pj]. As demonstrated in ref. |]Tl|, the 
resummation of threshold logs takes a particularly simple form, 

d^Tsi 



dwdxidz 



\Vub\ [Ksi {w, xi; as) f {z; as) + Dsi {w, xi, z; as)] , (16) 



7 



when one adopts the kinematical variables 



E 



w 



X 



Xi 



El 
'rriB' 



Ml 



:i7) 



The main point is that all IR logs are factorized in a function of the single variable z. Notice 
that z is only defined between 3/4 and unity. As already discussed, this formula can be 
generalized in the non-perturbative region corresponding to Mx ~ ExAqcd, where f{z) 
is related to the shape function defined in the effective theory by a short- distance factor. 
The 0{as) expressions for Ksi {w, xi; as) and D^i {w, xi, z; ag) can be found in ref. [0. The 
integration over w and xi involves only short-distance functions and can be easily carried 



out 25 1 . The resulting distribution reads 



dz 



iK^r r° [Ksi (as) f {z; as) + Dl^ {z- a,)] 



where 



si 



Gl ml m2_(m6) 



:i8) 



(19) 



1927r3 

We have chosen a somewhat unusual normalization factor T^i, in analogy with the radiative 



decay case. With our choice, the perturbative expansion of the total inclusive rate |2^ has 
a good convergence pattern. The coefficient is given by 



Ksi (as 



1 + 



a.Cf 



/403 _7r2\ +o/^2x ^_Qg^ 

Vi44 2 y ^ 



Let us now consider a change of the factorization variable of the kind 

1-z — >1- z' = 1- z + 0{l- z)^. 



(20) 



(21) 



It simply amounts to a rearrangement of higher-twist contributions and affects only the 
remainder function. Any ambiguity due to the choice of variable is related to the underlying 
leading-twist approximation and is common to all leading-twist treatments of the shape 
function^. The new factorized differential rate in terms of z' is 



dT 



si 



dz' 



Ksi{as)f{z', as) + D^iiz', as 



where 



T-,z' /■ I \ dz OLsCp 

E>si [z ; as) = —Dsi{z; as) + 



ln(l -z') + \ dz ln(l -z) + 



z' 



dz' 



z 



(22) 



(23) 



■^A first discussion of subieading twist effects in the sliape function region lias appeared recently ||2^. In 
B Xs-f they seem to be small for E-, < 2.3 GeV {x ^ 0.87). 
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Here we have left implicit the functional dependence of z = z{z'). The differential distribu- 
tions dVsi/dz and dVsi/dz' are two examples of observables which differ only by higher-twist 
effects. In principle, a comparison between observables of this type could provide some in- 
formation on the size of higher-twist contributions. 

As we are going to directly compare the semileptonic and radiative decays, it is preferable 
to choose the factorization variable for the two cases in a similar way. In i? — >■ Xs'^ the 
variable z = l- M'j^/{AE\) can be expressed in terms of the photon energy fraction x as 
z = 1 — Since we have factorized IR effects in the radiative decay in terms of a 

function f{x), it seems natural to adopt here a variable ^ such that 1 — z = (1 — ^)/ (2 — ^)^. 
Solving for ^ we then find 

i-e = ^, (24) 

where we have introduced the variable (0 < t < 1) 



t^^l_4(l_^) = ^l_Mi/Ei. (25) 

Unlike z, the variable ^ is defined between and 1. Using dz/d^ = ^/(2 — ^)'^ and the 
formulae given in [11|, it is straightforward to obtain Dsi{C,) = D^^i{^), which is also shown 
in fig. 3: 



Dsi (^; as 



IT 



r2 + 70 



(26) 



The higher-twist logarithms contained in the remainder functions are process-dependent 
and in general cannot be factorized in the function f{C,)- It is nonetheless interesting to 
note that the IR logarithm for ^ — 1 has a very small coefficient in eq. (^), much smaller 
than the corresponding coefficient of ln(l — z) in D^i{z). 
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One can also consider the partially integrated rate 

dr.. 



where at 0{as) we find 

■ 11725 - 6756 ^ - 5898 + 788 + 141 ^4 



(27) 



Esii^; as 



71 



+ 



5040 

^6-210^ + 135^2_5^3_g^4 
210 



ln(l - 



(28) 



If a cut on the hadronic mass Mx < is imposed, as often required in a realistic 
environment, the remainder function in eq. (^) is modified in the region 



m 



D 



nil 



The new expression is 



Dsl ii] Ois) |^<go 



3 C 



TT 



(1-0^ 



(7 + 6c) (e-2)e 



(29) 



70(2-O^ + 10c^ (1-0 + 7c (35 - 35^ + 8^^ 

35 



ln(l-0 



+ 



21-42^ + (21 + 114c2 + 16c3) -2c^ (57 + 8c) + 21c20 



+ 



(l-0' + c^e(2-3e + a + 



12(1-0' 

c^^ (22 -33^+ 15^2 _ 20)1 ln(l-0 



J (1-0= 



In the partially integrated case, the remainder function of eq. (pS]) becomes 
Esi\i<^o{i] as) = Esi\^>^o{^o; as) + A(^o; c) - A(^; c) 

for < .^0, with 



(30) 



A(e;c)= l^DsiiOde 
Jo 



s I 2 



TT 



(2-0 e 



56 c (1 -0 + 9 (22 - 22^ + 3 0) ln(l - 0' 



24(1-0^ 



7 _ 2 2c(14 - 28^ + 170 - 3 0) + 3 (33 - 66^ + 44 - 1 + 
4 ^ 6(1-0^ 
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(31) 
ln(l - 



360 c2 - 1) - 35 (96 - 96 ^ + 29 _ 14 c (246 - 246 ^ + 61 



315(1-0 



+ 2ci (e - 2) - + - ^) + ? -^H + se) _ ^) ^ ^ 

105(1-0^ 

More complicated experimental cuts can be applied taking as a starting point the triple 
differential distribution [lllf] and employing eq. (^31). 



4 Test of local duality and determination of \Vub/Vc 



cb 



In this section we construct a quantity that has a perturbative expansion in powers of as 
free from IR logs. In the shape-function region, this also guarantees the cancellation of 
non-perturbative effects related to the Fermi motion. Using the universality of the function 
/(x), it follows from eqs. (^) and (|18D that 



Ksi{as) 



dr 



si 



Krdiois) 



dV 



rd 



y*,\/,:pro, dx 



Drdix; as 



(32) 

Normalizing the experimental differential rates to the total semileptonic rate for B XJi', 
the previous equation can be rewritten as 



R = 



d jyo 
ljiXsI,u 



^^Dsii^;as 

9sl * 



^atird - — 



BR. 



9sl 



Drdix; as) 



Ksi{as 



77 



Krdias) 6q; 



ub 



(33) 



where BR^j^m and BR^rf are the branching ratios of the semileptonic and radiative decay, 
respectively. BR^^^c denotes the experimental value for the branching ratio of -B XJi^, 
BRsi^c ~ BRsi,tot = 0.1045 ± 0.0021 [^, while gsi = g{ml/ml) is the phase-space function 



for the semileptonic b ^ c decay, with g{x) = 1 — 8x + 8x^ — — 12a;^ Inx. Although it 
is sufficient to employ the tree-level expression for gsi in the above equations, there is an 
important uncertainty from quark masses. A very accurate determination of this phase- 
space factor, which could be adopted above, is gsi = 0.575 ± 0.017 The CKM matrix 
elements appearing in eq. (^) can be expressed in terms of the Wolfenstein parameters 
A = |Kts| — 0.22, p, and f/ as 



ub 



cb 



cb 



1 + A'(2p-1) + 0(A' 



0.97. 



(34) 



The r.h.s. of eq. ( P^D is a purely short-distance quantity calculable in perturbation theory 
and is independent of the point x = ^ where the l.h.s. is evaluated. Hence one expects to 
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find an experimental plateau in the range of ^ where the shape function provides a good 
description of the spectra. In practice, this allows a study of the extension of the shape- 
function region and a test of the validity of the local quark-hadron duality hypothesis. 
The numerical value of R depends on the CKM matrix elements. Let us introduce 



TT KJas 



1.94 ±0.16) X 10" 



(35) 



where we have estimated the theoretical error due to higher-order QCD corrections and to 
the top and charm masses on the basis of the results of the previous section, combining 
different uncertainties in quadrature. The ratio R is then given by 



R 



V,. 



ub 



5.15 X 10M1±0.08) 



V, 



ub 



VtbVtl 



(36) 



Of course, the theoretical error on R is not the only perturbative uncertainty affecting 
eq. (33), as higher-order corrections to the remainder functions Dgi and Drd might not be 



negligible, although they largely cancel against each other in the ratio R. In particular, we 
have seen that the remainder functions are relatively sizeable. The error on the r.h.s. of 
eq. (p3D induced in this way depends on the precise shape of the spectra and cannot easily 
be estimated. 

Provided a plateau can be experimentally identified, eq. (^) provides a stringent con- 
straint on the parameters of the unitarity triangle. In fact, it allows for a model-independent 
extraction of the ratio of CKM matrix elements l^b/Kfel- This is apparent from the second 
of eq. (p4D, which shows that | VjftVj*/V^f,p depends only weakly on p and f] (whose numerical 
values are taken from a recent global fit |2^). Therefore, to very good approximation, a 
measurement of the two spectra represents a measurement of |Ktfe/Kfe|- The master formula 
is 



ub 



V, 



cb 



BR 



sl,u 



dx 



\x=i 



C{as 



dx 



BR, 



rd 



(37) 



x=£. 



where C{as) has been defined in eq. (|35|) and 



6«BR 



sLc 



IT 



9si 



Drdii] ots) - ^Dsi{^; as 



K. 



si 



{31 



Whenever kinematical cuts have to be applied, the remainder functions have to be modified 
accordingly. In the case of a cut on the invariant hadronic mass, the relevant formulae are 
given in the previous section. 

We plot the function /i(^) in fig. 3. Quite remarkably, the two remainder functions cancel 
each other almost completely in the combination h{C,). In the range, ^ ~ 0.77, corresponding 
to the photon spectrum measured by CLEO |]3D|, the function /i(^) is small with respect 
to the measured spectrum. Higher order corrections to the remainder functions could 
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Figure 3: The functions h{^) and H{^) with and without a cut on the invariant hadronic mass 
(/if, = mb). For each case, the three lines correspond to values of Krd/ Kgi equal to our central 
value ±8%, see eq. i^. 



partially disrupt the cancellation only if they were very different in the matrix elements 
of the semileptonic and radiative decays. An additional source of uncertainty is the factor 
Krd/Ksi in eq. (|38|). The effect of the error on this factor given in eq. (^) is also shown in 
fig. 3. Although it is difficult to assess the impact of higher-order corrections on h{E,)i we 
can certainly neglect it in the final expression as a first approximation. These considerations 
are fairly independent of the presence of a cut on the hadronic invariant mass. Of course, 
it is straightforward to include in the analysis the small correction due to the fact that 
|Vtb\^*/V^b| is not exactly one. 

In a realistic analysis, it may be necessary to introduce some smearing over the spectra 
before applying eq. (^). Some smearing can be provided by partial integration over the 
spectra. In this case our ratio takes the form 



Vub 



BR. 



gsi 



EM) 



BRrd(x) 



BR. 



gsi 



Erdix) 



R, 



(39) 



where BRg^ and BKrd{x) are the partially integrated rates. The value of R is given in 
eq. (PBD. In analogy to the above discussion, one can vary the value of x = ^ where eq. (P^) 
is evaluated and look for the presence of a plateau. The master formula of eq. (|37D then 
becomes 



ub 



cb 



BR 



rd 



X 



BRrd(x) 



(40) 
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where again we have neglected 



6 a BR 



sLc 



9sl 



rd 



si 



(41) 



Indeed, as can be seen in fig. 3, as) is much smaller than the partially integrated 
branching fraction for B Xg'y (we recall that BRrd(0.77) 3 x 10""^ [0)- Although 
the uncertainty in the determination of H{E,) is large in relative terms, it seems unlikely 
that it adds much to the total uncertainty in the extraction of |V^b/V^fe|, which is therefore 
dominated by the error on C{as)- 

We stress that, since relation (|3^) is valid point by point, it is also possible to take into 
account a non-uniform experimental efficiency in a model independent way, provided that 
they are similar for the two rates. For example, in the ideal case of a common efficiency 
e(x = ^), eq. ( |40| ) will involve / dx' e(x') dSKrd,si I dx' . In the presence of different efficiencies 
for the semileptonic and the radiative decays, eq. (|40|) cannot be used in its present form and 
alternative startegies should be envisaged, using for instance the same equation translated 
in the A^-moment space or deconvoluting the efficiencies. 



5 Conclusions 

We have constructed a ratio of spectra of semileptonic and radiative B decays, from which 
non-perturbative Fermi-motion effects completely cancel out. This is a short-distance quan- 
tity with an as expansion not involving any IR logs. In spite of the semi-inclusive character 
of the spectra involved, the expansion is similar to those of inclusive quantities, such as the 
total semileptonic width. Combined with the experimental determination of the relevant 
spectra, our calculation allows for a test of the theory of the shape function and of the 
local quark-hadron duality hypothesis. If a kinematical "shape-function region" , where the 
shape function describes well the observed spectra, exists, the ratio R introduced in eq. (p^) 
should be independent of the variable x = ^. The quantity R involves the photon spectrum 
of radiative decays and the differential semileptonic charmless rate with respect to a certain 
function ^ of the hadronic mass over energy, Mx/ Ex- 

If indeed R is constant over a sufficiently extended kinematical region, our formalism 
permits a transparent extraction of |V^Mb/K:fe|. This is illustrated by a very simple formula, 
eq. (pTl), that links this CKM combination to the above two spectra. Thanks in part to a 
remarkable cancellation of NLO effects, eq. (|37|) is a very clean probe of the CKM structure. 
As this relation is valid point by point, the effects of experimental cuts and of a non- 
uniform experimental efficiency can in principle be incorporated. Assuming that the data 
confirm the hypothesis of small higher twist effects, a determination of \Vub/Vcb\ with a 
~ 5% theoretical error can be eventually obtained without using any model for the shape 
function. 
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The program we have just outhned can aheady be carried out using CLEO and LEP 
data, in parallel to the other |Ktfe/Kfe| analyses, based on the electron energy and the 
hadronic invariant mass spectra. Unfortunately, the available photon spectrum |^ is pro- 
vided in the lab frame; this effective smearing may hamper a precise comparison with the ^ 
distribution. Moreover, LEP has a small sample of events and a low signal to background 
ratio, and therefore the ^ distribution cannot be measured very precisely. We hope that 
some of these problems will be overcome at the B factories. Even with the present limi- 
tations, however, our method could provide a welcome check of the alternative techniques 
and we urge our experimental colleagues to implement it. 
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Appendix 



In this appendix we give all the relevant constants and functions necessary for the calcu- 
lation of the differential spectrum dTrd/dx. The non- vanishing functions dij{x) are given 
by 11,11 



dn 
di7 

d22 

d27 
d28 

dss 



x] 



x] 



X) 



x] 



X) 



1 

36 



d22 (x) 



1 

— ( 

6 



--d27 (x) , 



16c 



dt{\ - ct) 



/ dtRe 

9 Jo 

-^d27 (x) , 
1 

3 
2 

9 



1^12 (x) = -^d22 (x) , 
dl8 (x) = -^(^28 (x) , 

6 

G{t) 1 ' 
t ^2 ' 



[7 + X - 2a;' - 2 (1 + x) In (1 - X 

2 . /. , In (1 - X 
4 + x' + 4 (1 - x) — ^ 

X 



(42) 



27x 



-8 + 8x - x' - 2x^ + 2 (2 - 2x + x') In 



mr (1 — x) 



mi 
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The constants rj are given by 

1 



r2 , Re r2 = -4.092 - 12.78 (0.29 - mjrrih), 



6 

10 8 



9 ' 



Re rg = - (33 - 27r') . 



(43) 



The remaining rj with i = 3-6 have been calculated recently [Q, but their effect is very 
small because the corresponding Wilson coefficients are small. They can therefore be safely 
neglected. The constants ri and r2 depend very sensitively on the precise value of the charm 
mass. Following the arguments given in we use mc/mf, = 0.22 ±0.04 in ri 2 and obtain 
Re r2 = —4.99 ± 0.50. For completeness, we also give the value of 7^°^ entering eq. (p!OD, 



(0) / 208 416 176 152 6272 4624 32 32\ 



243 81 



243 81 



243 ' 3 ' 9 J 



(44) 



The functions ftj (x) entering the partially integrated spectrum read 



/ii (x) 
fir (x) 

fn {x) 

/28 {X) 

hi {x) 

/78 {x) 

/ss {x) 



36 



/22 (X) 



/l2 [X] = --/22 (x) , 
fl8 {x) = -^/28 (x) 

6 



— j(l-x)y^ dtil-ct) ^ + - dtil-ct) 



+ 7: 



_8c2 

"9" 

-lf27{x) 



rx/c 

X] I dtKe 
'0 



nl/c 

+ / dt{l-ct)Re 

J xlc 



Git) + - 



9 



;i - x) [31 + X - 2x^ - 3 (3 + x) In (1 - X 



(45) 



— [12 Li2 (x) - 12 (1 - x) In (1 - x) + 25 - 2ti'^ - 24x - x^] 



= — l-61n— (3-4x + x2 + 41nx) + 12 Li2 (x) + 24 lnx+ 
81 [ 

-3 (1 - x) (3 - x) In (1 - x) + (1 - x) (28 - 5x - 2x^) - 27r2|. 
The function G (t) in the integrand of /22 (x) and /27 (x) is given by 



G{t) 



-2 arctan^ ^' 4^t 



for t < A 



2W^i±|0_27rnn^*±|^-4 for t > 4. 



(46) 
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The above integrals can be performed analytically with the substitutions t = 4 sin^ y for 
t < 4 and t = A cosh^ y for t > A. The resulting expressions are quite long and are not 
reported here. For practical purposes, the compact integral representation given above is 
adequate. 



The relation between the dij and fij functions is simply given by 



dfij 

dx 



{x) . 



(47) 
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